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Abstract 



We compute the bulk limit of the correlation functions for the uniform 
measure on lozenge tilings of a hexagon. The limiting determinantal pro- 
cess is a translation invariant extension of the discrete sine process, which 
also describes the ergodic Gibbs measure of an appropriate slope. 



Introduction 

In this paper we study a well-known combinatorial model of determinantal point 
processes in dimension 1+1 (one spatial and one time variable). The model 
depends on 3 positive integers a, 6, c and is given by the uniform distribution on 
the finite set f2(a, &, c) of combinatorial objects that can be described in several 
equivalent ways: 3d Young diagrams (in other words plane partitions) in an 
a x b x c box; tilings of the hexagon of side lengths a, 6, c, a, 6, c by rhombi of 
three types; and a dimer model on the honeycomb lattice (see e.g. Cohn-Larsen- 
Propp |CLP| . Johansson [J], Johansson-Nordenstam JN]). Correspondences 
between different models can be found in the appendix. 

Following Johansson, we use yet another model where f2(a, 6, c) is identified 
with the ensemble of non-intersecting polygonal paths on the plane lattice. The 
latter model leads to determinantal point processes (put it otherwise, random 
point configurations) varying over time. The determinantal property, which 
is of crucial importance, means that dynamical (i.e. space-time) correlation 
functions of the model can be obtained as minors of some matrix which is called 
dynamical correlation kernel. The aim of the present paper is to compute the 
asymptotics of this kernel in the "bulk limit" regime. 

For a fixed time moment our model provides a random point configuration 
on the one-dimensional lattice — the so-called Hahn orthogonal polynomial 
ensemble. The dynamical model can be described as a chain of such ensembles 
with varying parameters, thus it can be called a dynamical Hahn ensemble. The 
dynamical correlation kernel is the so-called extended Hahn kernel; it has been 
previously obtained by Johansson [J] and Johansson-Nordenstam }JNI . The 
Eynard-Metha theorem [EM] provides the basis for computations of the kernel. 

The main result of our paper is computation of the asymptotics of the 
extended Hahn kernel in the bulk limit regime. The answer is given by the 
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translation-invariant kernel K(x, s;y,t) defined on Z 2 x Z 2 , for which a simple 
integral representation was derived: 

K(x, s;y,t) = & (f + cwf" 3 w^^dw. 

2iri J e - t 4, 

Here the integration is to be performed over the right side of the unit circle 
when s > t and over the left side otherwise, parameters <f> and c depend on the 
location in the bulk (see Theorem 1). 

The static version of the kernel (s = t) is the well-known discrete sine kernel 
onZxZ BOkOlJ. The obtained dynamical extension of the sine kernel is among 
the kernels constructed in |Bor] . The kernel K(x, s;y,t) is connected with the 
kernel obtained by Okounkov and Reshetihin in [OR] by a simple transformation 
(duality of particles and holes on the lattice). 

In the paper by Kenyon [Kenj local correlations in the model similar to the 
one considered by us in the similar limit regimes were studied. Using different 
techniques it was shown that the limit distribution is one of the so-called ergodic 
gibbs measures /i„ on the infinite hexagonal graph. The latter measures were 
classified by Sheffield in his thesis [She], and in the article [KOSj correlation 
kernels for the measures [i v were computed. Kenyon points out that his results 
should hold in the general case but in his paper |Ken| , full proof was given under 
the assumption that the so-called limit shape has no facets. In our case the limit 
shape does have facets, and it seems that our result does not follow from [Kenj . 

In the present article we do not employ the techniques of the papers by 
Kenyon, Okounkov and Sheffield and obtain Theorem 1 using a new method. 

Subtle analytical results about asymptotic properties of discrete orthogonal 
polynomial ensembles can be found in literature (see [BKMMJ, [j'J), but they 
deal with static rather than dynamical ensembles. The usage of the well-known 
technique of investigating asymptotics via a suitable integral representation of 
the kernel in the case of the extended Hahn kernel fails too. 

We solve the problem using the method first proposed in B02J. Although in 
|BQ2] again only static models were investigated, it turned out that the method 
can be adapted to the dynamical model of our interest. We use some relations 
between Hahn polynomials with varying parameters for the computation of the 
pre-limit kernel (similar approach was used in [BOlj ). 

The paper is organized as follows: In the first section, we define a stochastic 
process under study and compute its simplest characteristics. In Section 2, the 
determinantal form of the process is proved and its dynamical correlation kernel 
is expressed through the Hahn orthogonal polynomials. In Section 3, the main 
theorem is formulated and proved. Finally, in the fourth section we analyze 
obtained results and show their connection with theorems proved in [ULPj and 

The author would like to thank A. Borodin for proposing the problem, ad- 
vices and hinting the required form of relations for the Hahn polynomials. The 
author is grateful to G. Olshanski for many fruitful discussions and invaluable 
help in writing this article. 
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1 Basic definitions and properties 



1.1 Base model 

Consider the integer lattice 1? on the plane (i, x) and two sets of points on it: 
Xi = (0, i — 1) and Yi = (T, S + i — 1), where i ranges from 1 to some fixed 
N, while T and S are some fixed parameters (T > S). We are going to study 
collections of N non-intersecting paths passing through the lattice nodes, such 
that the i-th path connects Xi and Y{. Every path is a polygonal line each link 
of which is a segment with left end-point (t, x) <G Z 2 and right end-point (t+1, x) 
or (t + l,x + 1). In other words, the paths consist of horizontal segments and 
segments inclined at an angle of 45 degrees, thereby while moving along the 
path the ^-coordinate increases and the x-coordinate doesn't decrease. We will 
often call the ^-coordinate "time" . 

Introduce the probability space of all such families of paths and the proba- 
bility measure on it considering all families cquiprobable. 

Now let us consider an arbitrary time moment t and a set of points 
Zi,..., Zn with coordinates (t, z{) respectively (we assume that Zi are arranged 
in the ascending order). The probability of the event that our collection of 
paths passes through the points Zi is well defined. This probability induces 
a probability measure Pt{z\, . . . , Zn) on the TV-point sets of positive integers 
in a natural way. Next, the collections of paths under consideration define a 
stochastic process with discrete time t — 1, . . . , T, the state space consisting of 
iV-point sets of positive integers, and one-dimensional distributions P t . Let us 
denote this process by H t . Our aim is studying properties of this process. 

The above process is a Markov process: this is readily seen from the very 
definition of the Markov property (the past and the future are independent given 
the present). 

It is noteworthy that all paths ( and consequently all configurations of the 
process H t ) are contained inside a hexagon whose two sides are parallel to the 
t-axis, other two ones are parallel to the x-axis and the remaining two sides 
are inclined at an angle of 45 degrees. The corresponding illustration is given 
in figure 1 
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Figure 1. Collection of 6 non-intersecting paths. 

1.2 Elementary numerical characteristics 

We start with the explicit computation of the quantity P t (z\, . . . , z n ) 



Proposition 1. The number of collections of the non-intersecting paths con- 
necting points (ti,at) and (t2,&i), where i varies from 1 to N, is equal to 

Vbi — cij) 



detj jJ= i r 



,N 



Proof. This fact is well known. The trick which is used in the proof is described 
in [GY]. □ 

It follows from Proposition 1 and definitions that 
P t (zi, . . .,zn) 



dct ij= i i ... ! jv 


(z, + l-j)_ 


• detj^i,...,^ 


( T ~ t V 

KS+i-l-Zj) 




det lJ= i, 


..,N 


(s+i-j) 





Next we use Theorem 26 from |Kr| for the computation of the determinants. 
After some simplifications we obtain (in what follows we use the standard no- 
tation (a)i = a(a + 1) . . . (a + i — 1) ) : 



N 

n 



n ZjJ Li , l{t _ N _ ms _ N _ mT - t -s+ Zi )\ 

V (t + l) i _i(T-t + l) i _i(5-i + JV)!(T-S , + i-l)! ftl{T-t)\\ N 



n 

i=l 



(r + i)i_i(i-i)! 



T! 
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After further transformations the probability can be written in the form: 

Here Z is some normalization constant depending on T, N, S and t, while 
coordinate z[ is related to Zi simply through the shift by a quantity depending on 
the remaining parameters . All explicit formulas connecting different parameters 
in different representations of the weight function will be given in the next 
section. 
Thus, 

1 N 
P t (zi,...,z N ) = - Yl ( z 'i ~ Zj) 2 Y[ Wa 'P> M ( Z i)> 

l<i<j<N i=l 

where w a ^ t M( x ) — ^° + z!(m-z)! M ~ x ' x = 0, ■ ■ • , M, is the weight of Hahn or- 
thogonal polynomials (for the definition and properties of the Hahn polynomials 
see jKS] ) while parameters a, (3, M depend on T, N, S and t; as before, z- is 
related to zi by a shift of the origin. In what follows wt(x) denotes the weight 
function at the time moment t. 



2 Determinantal representation. The correla- 
tion kernel. 

2.1 One-dimensional distributions. 

Proposition 2. Consider an arbitrary non-negative weight w(x) on the set 
of non- negative integers. Assume that p$,pi, . . . are corresponding orthogonal 
polynomials. Further consider the probability measure on the N-tuples of integers 



1 N 

P t (Zx, . . . , Z N ) = — Y[ ( Z i~ Z j) 2 Y[ W i Z i)- 



l<i<j<N i=l 

Then the probability that the random N -tuple contains given distinct points X\ < 
x 2 < ■ ■ • < Xk is equal to 

det [K(xi,Xj)] 

where the kernel K is computed using formula 

N-i 



K(x, y) = y/w(x)w(y) 2J 



Pn(x)Pn{y) 
^ (Pn,Pn) 

Proof of this theorem can be found in (Ml 
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Such a probability distribution on the .ZV-tuples of integers is called the or- 
thogonal polynomial ensemble. Observe that one-dimensional distributions in 
our problem have precisely the same type as those in the last proposition. 

Now we return to the process H t under consideration. We are going to 
describe explicitly its one-dimensional distributions and the involved orthogonal 
polynomials and weight functions. 

In our model the weight function can be written in the form 



v ; Zi\{t -x + N - 1)\{S — x + N — 1)!(T - t - S + x)\ 

x'\{M-x')\ 

(here c is some constant that does not depend on x while x' differs from x by 
some shift). Note that the support of the weight function varies in time, i.e., 
when t increases by one either the support doesn't change or one point is being 
added to it from above or one point is being removed from below or both actions 
happen simultaneously (see Figure 1). Denote the support of wt by %t- 

Denote by Qk(x' , a, f3, M) the Hahn orthogonal polynomial of degree k with 
parameters a,j3,M. According to [KSj these polynomials are defined when 
x' = 0, . . . , M, k = 0, . . . , M, and they are orthogonal with respect to the 
weight 

(a + l) x >(/3+l) M - x > 
x'\{M-x')\ 

The Hahn polynomials can be expressed through the hypergeometric series: 

n ( i a »n ( -k,-x' ,k + a + (3 + l\ 

Q k {x ,a,/3,M) = 3 F 2 \ _ M a + l K 1 ) 

Denote by Ht(x) the Hahn polynomial Qk corresponding to the time moment 
t and shifted in such a way that its domain of definition coincides with Xt- 

There are 4 cases depending on the relations between t and other parameters. 
In each case a, 0, M and x 1 are expressed in terms of the parameters of our 
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model differently: 

(I) t<S + l, t<T-S+l: 

M = t + N -1, a = -S-N, f3 = S -T ~ N, x' = x, 

H?(x) = Q k (x>,a,0,M) = 3F 2 (- k >-*' A ^ +1 )(l) = 3 F 2 (:^^™)(1), 

0<x<M = t + N-l, < k < M = t + N - f; 

(II) S-l<t<T-S+l: 

M = S + N-1, a = -t-N, @ = t-N-T, x' = x, 

H?(x) = Q k (x',a,(3,M) = ,F 2 {- k ^ + a ^ +1 )(l) = M'-s-^-t-lX^ 

Q<x<M = s + N - 1, 0<k<M = s + N-l; 

(III) T-S-l<t<S + l: 

M = t + N-S-1, a = -T + t — N, /3 = -t-N, x' = T-t-s + x, 
H*(x) = Q k (x',a,f3,M) = ^{^^-^^^(l), 
-T + t + S <x' <t + N -1, Q<k<T + N -S-l; 

(IV) t>T-S-l, t>S-l: 

M = T-t + N-l, a = -T-N + S, = -S - N, x' = T-t-s + x, 
H k (x) = Q k (x',a,P,M) = 3F 2 (-^+ t +^; 2 ^-/ + + 1 1 )(l), 
-T + t + S <x' <S + N-1, 0<k<T + N-t-l. 

It is noteworthy that the expressions of polynomials H l k in terms of hyper- 
geometric series coincide in cases (I) and (II). The same is true for cases (III) 
and (IV). 

Let us introduce one more notation: 

Jn{ y/(HlW) ' 

Here (H„, H„) is the squared norm of the polynomial H l n with respect to 
the weight function wt(x). According to |K§] 

(Hi Hi) = (Q k (x', a, 0, M),Q k (x', a, 0, M)) 

= (-l) fc (fc + a + + l)jy+iQ3 + l) fc fc! 
(2k + a + 13 + l)(o + l) k (-N) k N\ ' 

The functions f^(x) form an orthonormal basis in the space l 2 (£<). 
Now the kernel K from Proposition 2 can be rewritten as: 

N-x 

K(x,y)=^2f n (x)f n (y). 

n=Q 

2.2 Dynamical determinantal property 

Let X(t) be a stochastic process taking values in the space of all subsets of some 
countable set. Let us say that X(t) is a determinantal process if there is a func- 
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tion K(x, s; y, t) such that for any n-tuple of distinct pairs (ari, ti), . . . , (x n , t n ) 
in the space-time 

Prob{£i e X(t 1 ),x 2 £ X(t 2 ), ■ ■ ■ ,x n £ X(t n )} = det [K{xuti;Xj,tj)\. 

i,j=l,...,n 

We will call the function K the dynamical correlation kernel. The above 
probabilities will be called the dynamical correlation functions. 

To prove that H t is determinantal we need the following abstract proposition. 

Proposition 3. Assume that for every time moment t we are given an or- 
thonormal system {/^} in linear space ^({0, 1, . . . , L}) (it may form the basis 
of the whole space or only of l^iLt), where L t C {0, ...... L}) and a set of 

numbers Cq, c* , . . . . Denote 

v t ,t+i{x,y) = y ^c t n f t n {x)f+ l (y). 

n>0 

Assume also that we are given a discrete time Markov process X t taking 
values in N -tuples of elements of the set {0, 1, . . . , L}, with one-dimensional 
distributions 

and transition probabilities 

det [vtj+ijxuyj)] det [/ffifaj)] 

N-l 

det[fU(^)] n 4 

n=0 

Then 

p n (xi,k 1 ;...]x n ,k n ) := Prob{xi £ X kll ...,x n £ X kn } 

= det [K(xi : hi] Xj, fej)] , 

i.j — l,...^n 

where 

N-l 

K(x,k;y,l)=J2 -Tkfi(x)f!(y),k>l; 

i=0 C i 

K(x,k;y,l) = -J2 C Y fiWfttv), k < h 

i>N 

= 1, cW = cf • c^ 

This proposition is just a slight extension of the theorem proved by Einard 
and Metha in the paper [EMJ . In their article this proposition was proved under 
the assumption that the system {/^} forms a basis of the whole space. But this 
assumption is not essential as we always can extend an orthonormal system to 
a basis of the whole space and then set the coefficients c* corresponding to the 
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extra vectors equal to zeros. This operation doesn't change neither transition 
probabilities nor one-dimensional distributions of the process. 

Another similar formulation of the proposition was used in [BOlj . 

It is noteworthy that the dynamical correlation kernel is not canonically 
defined. Only dynamical correlation functions are invariants of the process and 
they don't change when we transform the correlation kernel in the following 
way: 

K*(x, S ;y,t) = K(x,s;y,t)-^^-, 

here F is arbitrary function. It is evident that such transformation doesn't 
change the values of determinants. 

Lemma 4. One-step transition probabilities of the process H t can be computed 
by formula 

Pt,t+i(xi, ■ ■ ■ ,x N ]yi, ■ ■ -,Vn) ■= 

= Prob{H t+ i = (z/i, . . . ,y N ) | H t = (x u . . .,x N )} = 

mvj-vi) n (N+s-xi-i) n fa+T-t-s) 

i<j i: yi—Xi-\-l i: yi—Xi 

(t - t) N ■ n (xj - Xi) 

provided that each difference yi — Xi is equal to zero or to one, otherwise the 
probabilities are equal to zero. 

This formula can be checked by direct computation starting from the def- 
inition and using Proposition 1 together with the determinantal formula from 
|Krj used in Section 1.2. 

Proposition 5. The process H t meets the assumptions of Proposition 3 and 
thus it is determinantal. 

Here the orthogonal system from the formulation of Proposition 3 is the 
system of functions /„ defined in Section 2.1 while coefficients c\ are given by 
the following formula 



C* = 4 / I 1 - ~ " 7T III- 



t + N J V T + N-t-1 
Proof. We will need 3 lemmas. 

Lemma 6. The following relations hold for the Hahn orthogonal polynomials: 

(*) xQ k (x-l;a,[3,M-l) + (N-x)Q k (x;a,f3,M-r) = MQ k (x;a,f3,M), 
(**) xQ k (x- l;a+ 1,(3- 1,M) + (-x-a - l)Q k {x; a + 1,0 - 1, M) = 

= -(a+l)Q k (x;a,0,M). 
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They are equivalent to the following relation for the function 3F2: 

("cr + ? y>-^ (t;;T) w - ( _ \r) «■ 

Proof. The following formula holds (see |KSj ): 

nt a m\ t? (~ k > ~ x ' k + a + ^ + 1 V-n 
Q fe (^a I /3 ) M)= 3 F 2 ^ _ M;a+1 J(l) 

^ (-M)i(a + l)(i! 

It implies 

xQ k {x - 1; a, [3, M - 1) + (AT - a;)g fc (x; a, /3, M - 1) - iVQ fc fo a, (5, M) = 



k 



y, (-k)i{k + a + f3+ 1), 



=0 

A' 



(-x + l)i (-x)i (-x) l 



x-r + (M- x)-^ ^— - M , 

(-M + l)i V '(-M + l), {-M)i 



i=0 



The second relation is readily seen to be equivalent to the one already proved. 

□ 



Lemma 7. 

M 

i--i.)'-ia- f ' ./ riw.-i! ;-m );,/,! 



(2fc + q + /? + l)(a + l) fc (-M) fc M! ^ 



(a + l) x (/3 + l) M -x 
jc!(M-a:)! 

This formula is called "dual orthogonality relation" and can be found in 
El- 
Lemma 8. 



Iw t +i(y) 



5 + T - 1 - g + x xy 

f°x + l ~1 Jr. , 7\ t" \ / m , at , TV a:' 



,/(i + jv)(r + Jv-t - 1) + ,/(t + Jv)(r + Jv-t-i) 

/fere coefficients c\ and functions ft(x) were defined in the formulation of Propo- 
sition 5, Wt(x) is the weight function of the Hahn orthogonal polynomials. 
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Proof. To prove this it suffices to recall the definition of fl(x), substitute the 
value of the weight function and all parameters, then in cases (I) and (11) express 
ff, +1 (y) in terms of fl(y) and fl(y — 1) using the relation (*), and finally use 
the last lemma. In the cases (III) and (IV) one should express fl(x) in terms 
of fl +1 (x) and fjf +1 \x + 1) using (**) and again use the last lemma. □ 

Let us return to the proof of the proposition. 

The fact that the one-dimensional distributions can be expressed in the form 
that is required for the application of Proposition 5 follows from Proposition 2, 
where one takes k — N (actually, in this case Proposition 2 is quite obvious) 
and the following observation: 

(. . det N [fUix^f = d„ v [fUxj)] ■ . . = det [fUM] 

= det([fU( Xj )] ■ [fU&i)}) = . . det \K{ Xj , Xi )] 

(here K(x,y) is the kernel defined in Propostition 2). 

It remains to show that the transition probabilities also have the required 
form, i.e: 

p , x det [v tit+1 ( Xi , yj )]det [f-^(y 3 )] 
Pt,t+i(xi, ...,x N ;y 1 ,...,y N ) = — 

det [ftM)] n ci 

n=0 

(all matrices under the sign of determinant have order TV x N and we will not 
indicate this in the sequel). Here 

v t , t+1 (x,y) = ^clfl(x)f t k + \ y ). 

k>0 

By virtue of already proved facts we obtain: 
v t , t +i{x,y) = 



I w t {x) 



S + N-l-x ru T-t-S + x 
— fry _| 

y/(t + N)(T + N-t-l) x+l y/(t + N)(T + N -t-l) x 



det [ft+l{ yj )\ = v / P t+1 (y 1 ,... 1 y N ), 



det [fi-^Xjj] = y/P t ( Xl ,...,x N ). 

Substituting the expressions of the weight function and of the probability of the 
A^-tuple we get: 

det [v t ,t+i(xi, yj)} det [f^iyj)] 

N-l 

det[/t!fe)] II 4 

i=0 
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det [(S + N-x- l)8% +1 + (T-t-S + x)8 v x \] „ Vj 



(T-t) N f*j Xj - x % 

Let us assume that and yi are arranged in the increasing order. 
There are two cases: either every difference y, — a;, is equal or 1 or this is 
not true. 

In the latter case the matrix under the sign of determinant has a zero row 
or column. Consequently, the determinant is equal to zero. On the other hand, 
the transition probability is equal to zero, too. 

In the former case the matrix under the sign of determinant is block-diagonal. 
Moreover, its zeros are arranged in such a way that every block is either an upper 
triangular or lower triangular matrix. We may conclude that the determinant 
of the matrix is equal to the product of its diagonal elements, and the value 
of (* * *) is equal to the transition probability of the process H t computed in 
Lemma 4. 

Hence the process meets the conditions of Proposition 3 and we conclude 
that the process is determinantal. □ 



3 Passing to the limit in the correlation kernels: 
the operator method 

3.1 Statement of the result 

We are interested in the limit distribution (limit correlation kernel) in the fol- 
lowing limit regime: 

Let us fix numbers S, T, N, t, x. Assume that auxiliary parameter p — > oo 
and 

S = pS + o(p), T = pf + o(p), N = pN + o(p), t = pt + i, 
s = pt + s, x = px + x, y = px + y. 

This limit has a simple geometrical interpretation: we fix the proportions 
between the parameters T, S, N of the hexagon containing the paths and a 
point (t, x) inside it. Further we enlarge the hexagon by means of homothety 
with center at (i, x) and coefficient equal to p. Thus, the number of the paths 
increases, and we are focusing on the local picture near the point (t, x). 

Recall the formula for the pre-limit kernel: 

K(x,s;y,t)= - £ [f[ A f!{x)ft{y), s <t 

i>N \j=s ) 

N t' [ft y\ fl(x)ft(y), s>t. 

i=0 \ j=t c » / 
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The functions ff(x) and the coefficients c\ were defined in Section 2. 

Theorem 1. Let all the parameters of the original model (S, T, N, t, s,x,y) tend 
to infinity in such a way that 

S = pS + o(p), T = pf + o{p), N = pN + o(p), t = pt + i, 
s = pt + s, x = px + x, y — px + y, 

where the auxiliary parameter p — > oo. Then the limit point process on Z 2 exists. 
It is determinantal and has translation-invariant kernel 

K(x, s; y, i) = -L I (1 + cw)^' s w^'^dw. 
2m J e - i<t , 

Here the integration is to be performed over the right side of the unit circle when 
s >t and over the left side otherwise, 




x(S + N-5:) 
(T -t- S + x)(t + N - x)' 



and angle <j) is given by formula: 

-N(N + f) + (-x + S + N)(t + N-x)+x(f + x-S~ t) 
(p — arccos . 

2yjx(-x + S + N)(t + N — x)(x + f-S-t) 

If the expression under the arc cosine sign is greater than 1, then we set = 0. 
// the expression is less than —1, then <j) = 

Existence of the process given by the limit correlation kernels follows easily 
from the convergence of these functions. The proof of this fact can be found for 
instance in |Borj (Lemma 4.1). 

The proof of the convergence of the correlation kernels will be made in two 
steps. 

3.2 Statical case 

First assume that s — i — 0, i.e. we are investigating only the "statical" distri- 
bution in some fixed time moment. 

We use the method first proposed by A. Borodin and G. Olshanski in the 
paper [BQ2] , 

We extend the functions fj{x) belonging to the space hiX-t) to the functions 
defined on Z by setting them equal to outside %t ■ Now we may regard all the 
functions /' as belonging to one and the same space fo(Z). 

We have 

N-l 

K t (x,y)=^2f i (x)fl(y). 

i=0 
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Observe that K t (x, y) is the matrix element of the projection operator onto the 
subspace spanned by the first N orthogonal functions /*, i = 0, . . . , iV — 1, in 
the space Z 2 (Z). 

It turns out that finding the limit of the operator is easier than computing 
the limit of the matrix elements. Note that functions f„(x) are eigenvectors 
of some difference operator (it will be explicitly given below). The projection 
operator can be regarded as the spectral projection on the segment containing 
the first N eigenvalues of the difference operator under investigation. Now, to 
find the limit of the spectral projection operators we will take the limit of the 
difference operators. Note that both the difference operator and the spectral 
segment are varying simultaneously. 

To justify the limit transition we use some theorems from functional analysis. 

Consider the set l\(7£) of the finite vectors from £2(2) (i.e. algebraic span 
of the basis elements S x ) as a common essential domain of all considered differ- 
ence operators. It will be clear from the following that the difference operators 
strongly converge on this domain. It follows that operators converge in the 
strong resolvent sense (see |RSj . Theorem VIII. 25). The last fact, continuity 
of the spectrum of the limit operator and Theorem VIII. 24 from RS] imply 
that the spectral projections associated with the difference operators strongly 
converge on the set of finite vectors to the limit spectral projection associated 
with the limit difference operator. 

Now we will carry the explicit computation through. 

The following difference relation holds for the Hahn polynomials (see |KSj ): 




k(k + a + j3+ l)Q k (x';a,l3,M) = (x' + a + 1)(M - x)Q k (x' - l;a,/?,M)+ 
+ {(x' + a + 1){M -x) + x'(x' -p-M- l))Q k (x'; a, (3, M)+ 



Let us rewrite this relation in the terms of orthonormal functions fl(x'). It 
must be emphasized that here the variable x 1 is used instead of the variable x 
(It was shown in Section 2.1 that in cases (I) and (II) x = x', while in cases 
(III) and (IV) x^x'). 

fc(fc+q+/3+l) ft[ x 'j _ y/{x' + l)(.M-x')(l3+M-x')(.a+x') j. t + ^_ 

(x'+a + l)(x'-M)+'x(x'-p-M-l) ft( x >} _|_ ^x'(M-x' + l)(l3+M-x'){a+x') _ ^ 

We divided the relations by p 2 to ensure the existence of the limits that will be 
considered below. 

Denote by H the difference operator given by the right side of the last equa- 
tion for those x for which the Hahn polynomials are defined and, for instance, 
as identity operator for other x. The matrix of the operator H is given by 




+ x'(x' -/3-M- l)Q k (x' + I; a, /?, M). 




K t (x,y), x,yeX t 
51*. for other x,y. 
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It is readily seen that H is a self-adjoint operator. 

We are interested in the projection on the eigenvectors corresponding to the 
part of the spectrum for k G [0, . . . , N— 1], i.e. corresponding to the eigenvalues 

fc(fc + Q + /3 + 1) . ke[0...N-i]. 
p 

Observe that the full spectrum of the operator H consists of the negative num- 
bers that can be obtained by last formula when k ranges from to M and the 
value 1 that occurred because of the formal continuation of the operator. In 
our model a and (3 are rather great in magnitude negative numbers and the 
function k(k + a + (3 + 1) is negative and monotonic for the permissible values 
of k. Consequently, we can replace the projection on the eigenspaces associated 
with the discrete set of points of the spectrum by the projection associated with 
some segment containing them. Thus the desired projection P is the spectral 
projection associated with the following segment: 

■ (N~l)(N + a + l3) n 
P 2 

Substitute the values of a, (3 and other parameters and express H and the 
projection P through p and the parameters S, T, N, t, x: 

S = pS + o(p), T = pf + o{ P ), N = P N + o(p), 
t = pt + t, x — pi, + x. 

We obtain the family of the operators P p in the space ^(Z) where x is the 
variable. Our aim is finding the limit 

P = lim P„. 

p — >oo 

For the application of Theorem VII. 24 from |RSj it is necessary that the 
family of the difference operators H p obtained by substituting the parameters 
of the limit regime converges to some limit H in the sense of the strong operator 
convergence on the space of the finite vectors. This convergence follows from 
the convergence of the matrix elements because in our case all operators H p are 
given by tridiagonal matrices. 

Taking the limit of the matrix elements we derive that H is a symmetric 
operator which is given by tridiagonal matrix that has the number 

A ,. (x' + a + l)(x'-M) + x'(x'-l3-M-l) 
A = hm '— = 

p^oo p z 

= -(§ + N- x)(t + N - x) - x{x + f - S - t) 
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on the diagonal and the number 



B = lim 

p — >oo 



y/(x' + 1)(M - x'){j3 + M - x'){a + x') 



= lim 

p — >oo 



y/{x'){M - x' + l)(/3 + M — x')(Q! + x') 



(S + N — x)(t + N — x)x(x + f-S-t) 



above the diagonal and below it. 

It is noteworthy that although the parameters a, /3, M and x') were ex- 
pressed through the characteristics of the model in the different ways (we had 
4 cases) here these differences disappear and the limit difference operator is the 
same in all cases. 

The limit spectral segment can be easily computed also: 



lim 

p — >oo 



(N-l)(N + a + 0) 



,0 



-N(N + T),0 



Denote by I\ a M the indicator of the segment [a, b] 
Observe that the following identity holds: 



P = I\ 



-JV(JV+f'),0] (^0 _ I\ -N(N + t)-A 



H - A 
2B 



H 2g A is the operator which is given by the matrix with the number ^ above 
and below the diagonal and zeros in all other entries. 

To compute explicitly P we make the Fourier transform ^(Z) — > L2(S 1 ) 
where S 1 is the unit circle in C. 



The Fourier transform of 



H-A 
2B 



is the operator of multiplication by the func- 



tion = Observe that the last operator has purely continuous spectrum 
that fills the segment [—1,1]. Hence, the spectrum of the operator H is purely 
continuous too (this condition is necessary for applying Theorem VII. 24 from 

USD- 
Next, we substitute A and B. The coordinate of the left end-point of the 
spectral segment is equal to 

~N{N + f) + (-x + S + N)(t + N- x) + X(f + x-S-i) 



2^x(-x + S + N)(t + N — x)(x + T-S-t) 

while the coordinate of the right one is equal to 

(-x + S + N){t + N -x) + x(f + x-S- t) 
2Ji{-x + S + N)(t + N — x){x + f-S-t) 
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Note that the value of the last expression is greater than one. Therefore, 
the Fourier transform of the spectral projection P becomes the operator of 
multiplication by the characteristic function of the right arc of the unit circle 
contained between the angles — and <f> where 

-N(N + f) + (-x + S + N)(t + N-x)+x(f + x-S- t) 
= arccos — 



2yjx{~x + S + N){t + N- S-yjc + f-S-t) 



Finally we perform the inverse Fourier transform and find the matrix of the 
operator P which coincides with the desired correlation kernel 

i r % * i r 1 * _-_ 

KAx,y) = K(x,y) = — - f w v x 1 dw = — - & w x y 1 dw 
ziri J e -i<t> 2m J e -i<f 

(the integration is to be performed over the right side of the unit circle) . 

The last integral can be computed explicitly and it coincides with the discrete 
sine kernel: 

K(x,y) = Sin f X -f . 

ir(x - y) 

Observe that applying the same method we could have found the limit of 
the kernel 

k t {x,y) = -Y,fKx)fi{y)- 

i>N 

This limit is equal to 

K i {x,y)=K{x,y) = —j ^w x - y - l dw, 

where the integration is to be performed over the left side of the unit circle. In 
the space ^(S 1 ) we obtain the operator of multiplication by the characteristic 
function of the left arc of the unit circle. 

3.3 Passage to the dynamical kernel 

Now assume that s < t. We have 

K(x,s;y,t) = -J2 f f[ 4 ) KWfHv)- 

i>N \j=8 J 



Fix s and t and let Q be the operator in li^E) which is given by matrix 



K(x,s;y,t), xeX s ,yeX t , 
for other x, y. 



This operator can be expressed in the following way: 

Q = Vt- U t -i ■ U t -2 ■■■u s . 
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Here Vt is the operator corresponding to the kernel K t (x,y) considered in the 
end of the last section, this operator is given by matrix 



- E f!(x)f!(y), x,yeX t , 

i>N 

for other x, y, 



and Uh is the operator mapping every function f^(x) into c' i / 4 /l+1 (x) (if in the 
time moment h + 1 one less polynomial exist, then the function corresponding to 
the polynomial of the highest power is mapped to zero; this definition agrees with 
the fact that corresponding coefficient c' 1 is equal to zero). On the orthogonal 
complement to the span of the functions //* set Uh equal to zero. This operator 
can be given by the matrix 

(E4fHx)f- +1 (y), xeX h ,yeX h+1 , 
U h (x,y) = < *>o 

I for other x, y. 

Now we express all parameters in the terms of characteristics of the limit 
regime and obtain the family of the operators depending on the parameter p. To 
compute the limit Q p when p^oowe take the limit of every factor separately 
and then use the following argument: all factors strongly converge on the dense 
set of the finite vectors, the norm of each factor is not greater than 1 (For the 
operators Uh it follows from the fact that < < 1), consequently, operators 
strongly converge on the whole space liil^). The operation of multiplication is 
continuous in the strong operator topology if norms of all factors are uniformly 
bounded. This argument implies the strong convergence of the operators Q p . 

Observe that the matrix Uh have already been computed in the proof of the 
determinantal property of the process: 

Uh(x,y) = v h ,h+i(x,y) = 

S + N — 1 — x cV T — h — S + x 




y/(h + N)(T + N-h-l) x+ y/(h + N)(T + N-h-l) * 



(S + N-x-l)(x + l) ( gy l(T-h-S + x)(h + N -x) 



(h + N)(T + N-h - 1) x y {h + N)(T + N -h-1) 

(x eX h ,y e X h+1 ). 

Next, we substitute in the last expression all parameters of the limit regime 

S = pS + o(p), T = pf + o(p), N = pN + o(p), 
t = pt + i, t = pi + s, x = pi + x 

and obtain the operators Uh, P - Taking into account that s < h < t we send p 
to infinity. We obtain the operator which is given by the matix 



u { x,y) = 4 J S s t N r^K + 4J( T -r s - + - £)ii + N -- £) . 

X+1 V {t + N)(T + N -t) X V {t + N){T + N-t) 
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Now we make the Fourier transform and multiply all factors. The result is 
the limit operator Q.It is the operator of multiplication by the product of the 
function 



l(T-t-S + x)(t + N-x) / (S + N-x)x 



\]j (t + N)(T + N-t) V (t + N)(T + N -t) J 

and the indicator of the left side of the circle contained between the angles (f> 
and 2ir — <f>. 

Finally, we perform the inverse Fourier transform and derive the desired 
correlation kernel: 

K(x, s; y, i) = 




(T-t-S + x)(t + N-x) 1 / {S + N-x)x 



{t + N)(f + N-t) z V (t + N)(f + N - i)^ 

■ z y-^dz w = /z 




(T-t-S + x){t + N-x) / (S + N-x)x 



{t + N){T + N-t) V (t + N)(T + N -t) / 



t-s 



(T — t — S + x)(t + N — x) \ 1 
{t + N)(f + N-t) ) 



f (l + cw)w x v 1 dw. 



The factor outside the integral can be omitted as it corresponds to a conjuga- 
tion of the kernel which doesn't change the correlation functions. We obtained 
the desired expression from the formulation of the theorem. 

In the case s > t reasoning is similar and we will omit some details. Consider 
the operator Q which is given by the matrix 



G(x,y) 



K(x,s;y,t), xeX s ,yeX t , 
for other x, y. 



It can be expressed in the following way: 

g = v s - 1 ---v t+1 -v t -n. 

Here Vh is the operator defined on the first N functions //* by the formula 
Vh(fi) = • As above we extended it to the operator on the whole space. 
Denote by 1Z the orthogonal projection onto the span of the first N functions 
ft- 

Next we express the operators through the parameters of the limit regime 
and obtain three families of the operators: Q p , Vh, P and 1Z P . 
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Using the definition of c\ we conclude again that norms of the operators Vh, p 
and 1Z P arc uniformly bounded. 

TZ, which is the limit of the operators 1Z P when p — > oo, has been already 
computed in the previous section. This operator is equal to the operator of the 
restriction to the right arc of the unit circle contained between the angles —<f> 
and (j). 

For further reasoning we need the following lemma: 

Lemma 9. Consider two families A p and B p of the operators in the Hilbert 
space. Assume that when p — > oo norms of all operators are bounded by some 
constant, A p — > A and B p A p — > P in the sense of the strong operator conver- 
gence, P is an arbitrary orthoprojection and E is its image. Also we assume 
that E is invariant space for the operator A and a restriction A E of the last 
operator is invertible. Then on the subspace E operators B p strongly converge 
toAz 1 . 

Proof. Consider an arbitrary vector x belonging to the E. It follows that one 
can find such y that x = Ay. Then we have: 

B p A p y -> y, 

A p y -> x, 

B p x = B p (A p y + {x- A p y)) = B p A p y + B p (x - A p y) ->y + = y. □ 

Let us apply the proved lemma taking A p = Uh+i, p and B p = Vh, P - The 
product B p A p is the operator which is given by the matrix that coincides with 
the kernel studied in the statical case and, thus, converges to the projection It. 
It was shown above that the operators A p converge to the operator U which is 
the operator of multiplication by the function 

j (f-t-S + x){t + N-x) | / (S + N - x)x 

V (t + N)(f + N-t) \j (t + N)(f + N-t) 

in the L2(5 1 )-realization. 

The restriction of this function to the unit circle is non-zero and invertible 
when z ^ — 1 while it may be equal to zero when z = — 1. But the limit 
projection TZ is the operator of multiplication by the characteristic function of 
some arc which doesn't contain the point z = — 1. Consequently, restriction of 
the operator U to the image of projection TZ is invertible and we may apply the 
last lemma. 

We obtain that (in L2(-S' 1 )-realization) limit of the operators Q p multiplies 
vectors from the image of the limit projection TZ by some function. This function 
is inverse to the function which the operator U is multiplying by. On the other 
hand, 1Z P — > for vectors from the orthogonal complement to the image of 
the projection 1Z, consequently, as norms of all operators Vh, P are bounded, so 

Q P ^0. 

Further reasoning is identical with the case s < t. Theorem is proved. 
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4 Analysis of the results 



4.1 Frozen regions 

Now we will take a look more closely at the obtained one-dimensional correlation 
function which can also be called the density function. 

Recall that all points of the collections of non-intersecting paths under con- 
sideration form a hexagon whose sides are parallel to one axis, other two ones 
are parallel to another axis and third pair of the sides is inclined 45 degrees with 
respect to the axis. 

In the paper [ULPj tilings of the hexagon by rhombi were studied. These 
tilings can also be interpreted as collections of non-intersecting paths. It was 
proved that so-called "frozen" (or "arctic") regions exist. In these regions the 
limit distribution is non-random and tiling consists of rhombi of only one type 
(for the exact wording see Theorem 1.1 in [CLP]). We will show below that 
our results agree with this theorem and the shape or the frozen regions can be 
obtained also from our dynamical kernel as regions of such values of the parame- 
ters of the limit regime that local fluctuations disappear and the limit process is 
trivial (either all points of the lattice belong to the random configuration with 
probability one or no points at all belong to the random configuration with 
probability one). 

Let us fix the parameters of the hexagon N, T, S in our model and study 
the dependence of the density function on the point (t,x). 

Proposition 10. The limit distribution is non-trivial only inside an ellipse, 
which is tangent to the hexagon defined by the parameters N, T, S while outside 
this ellipse the density function is equal to zero or to one. 

Proof. Let us compute K{x,t;x,t) which coincides with the density function. 



For some values of the parameters either 4> = or <fi — n, it means that the 
density is equal to or 1. In our terms the desired case is equivalent to the fact 
that in the expression defining the angle (/>((/> = arccosD) the parameter D is 
not less than 1 in magnitude. We can rewrite it as D 2 > 1. After simplification 
last inequality transforms into 



f 2 x 2 + {S + N) 2 t 2 + 2xt{NT -Sf- 2SN) + 

+ 2t{SN 2 - NTS - N 2 f + S 2 N) + 2i(NTS - NT 2 ) + 

+ N 2 (f-S) 2 < 0. 



Interior of some ellipse is given by this expression. One can easily check that 
this ellipse is tangent to the borders of the hexagon under consideration defined 
by the parameters N, T, S. 

In the paper CLPJ the frozen region was the same, it was the interior of the 



We have: 




tangent ellipse. 



□ 
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4.2 Connection with Okounkov-Reshetikhin model 

In this section we will show the connection between our results and results of 
the paper [OR] . In the last paper some point process was studied. It had the 
correlation kernel equal to 



here z is some complex number depending on the parameters of the model. 
This number is less than 1 in magnitude and has positive real component. The 
integration is to be performed over the right arc of the circle when t > s, 
otherwise, over the left arc. 

Proposition 11. One can choose a linear transformation of the lattice Z 2 in 
such a way that processes given by our kernel and by the kernel from \OR\j are 
connected by the particles-hole involution, i.e. random configurations given by 
the former kernel are complementary to the ones described by the latter kernel. 

Proof. Observe that choosing five parameters of our model we can obtain arbi- 
trary values of the kernel parameters < c < 1 and < </> < it. 

Now take a look at the kernel from [OR] . Let us change the variables in the 
integral to transform the path of integration into the arc of the unit circle (note 
that arising factor outside the integral can be omitted as this factor doesn't 
change correlation functions). Next, we introduce the linear transformation of 
the lattice 

s' - t' = s - t, 

x' - y' = x - y + (s- t)/2. 

Indeed this transformation maps the lattice considered by Okounkov and 
Reshetikhin into Z 2 . The kernel becomes 



(the integration is to be performed over the unit circle, the choice of arc is made 
in the same way as above, < c < 1) . For s ^ t we make a change of variables 
in the integral z = —w. 



Here the integration is performed over the right arc of the unit circle when t < s 
and over the left one otherwise. The factor (—l) x ~ v can be omitted again. This 
kernel differs from the kernel of our model by the sign for s ^ t. While for 
s — t the kernel differs from the desired one by the arc of the integration. To 
obtain the required arc we should subtract the integral from the residue of the 
integrand at the point 0, i.e. subtract from 1 for x = y and subtract from for 





(1 — cw) 
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x 7^ y. Thus, two kernels are related in the following way: K\ = 8 x=ytS=t — K%. 
But it is exactly the transformation of the kernel in the particle-holes involution, 
(for instance, it was proved in Appendix A. 3 in [BOkOlj ). 

Note that when the parameter of the integral c > 1 in our model we can 
change the variables: w = 1/z. The constant c becomes 1/c and also index 
of power of z changes from x — y — 1 into x — y + (s — t) — 1. However, this 
transform is only another linear transformation of the lattice. Consequently, 
the proposition is proved in that case too. □ 

This result could have been foreseen, because the models are closely related. 
In [QRj the ensemble of 3d Young diagrams without any restrictions was studied 
and one of the interpretation of our model is the ensemble of 3d Young diagrams 
contained in the box. Other argument in favor of this result is hypothetical 
universality of the obtained kernel stated in |OR| . One of the additional proofs 
of the universality was derivation of the same kernel in [Borj as bulk scaling 
limit of the cylindric partitions. 

5 Appendix. Possible interpretations of the 
model under consideration 

The model under consideration has lots of different combinatorial interpreta- 
tions. We will show some of them. 

Consider a collection of non-intersecting paths we studied through the paper. 



Figure 2. 3 non-intersecting paths. 

Now we incline at an angle of 30 degrees with respect to the horizontal axis 
segments of the polygonal lines which were inclined at an angle of 45 degrees. 
And replace the horizontal segments by the segments inclined at an angle of 30 
degrees to the axis again, but this time downward directed. This new picture 
can be interpreted as a family of the paths on the surface of 3d Young diagram 
in the box: 
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Figure 3. Non-intersecting paths on the surface of 3d Young diagram. 

The last image can also be regarded as a tiling of the hexagon by rhombi of 
3 types. 
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